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CHAPTER  I 


INTRODUCTION 


One  of  the  basic  problems  in  quantum  theory  is  the  solution 
of  the  time-independent  SchrBdinger  equation 

where  the  Hamiltonian  ^  is  bounded  below  and  has  associated  with  it 
a  set  of  solutions  (E,)  corresponding  to  the  discrete  levels  lying 
below  the  continuum. 

This  equation  cannot  be  solved  exactly  except  for  a  few  cases 
corresponding  to  very  simple  systems.   Instead  the  problem  is  solved 
by  adopting  a  suitable  method  for  obtaining  an  approximate  solution. 

One  particularly  useful  method  for  treating  this  problem  is  the  parti- 

1-6 
tioning  technique.     This  approach  has  the  interesting  feature  that 

it  contains  many  of  the  more  conventional  methods  such  as  the  variational 
principle  and  perturbation  theory  as  special  cases. 

Using  the  partitioning  technique  one  may  construct  a  function 
C,  =  f(C)>  where  £  is  a  real  variable,  having  real  solutions  €■,-,» 
Sl2*^lV"''^l     ^^^^  ^^^   property  that  each  pair  {£,£,.)     will 
bracket  at  least  one  eigenvalue  E  of  the  Hamiltonian  Ji.  .   This  means 
that  if  £    is  chosen  as  an  upper  bound  to  an  eigenvalue  E,  at  least 
one  of  the  solutions  g  ..   will  be  a  lower  bound  to  E.   The  problem  of 
constructing  this  function  and  a  way  of  determining  which  of  the 


solutions  €  ■,  .    is  the  lower  bound  to  E  will  be  studied.   The  lower 
bound  to  E  is  taken  to  be  the  largest  of  the  solutions  Ci-'  (^~^'^» 
...,k,k^g),  providing  a  lower  bound  to  E.   Furthermore,  it  will  be 
shown  that  lower  bounds  to  the  solutions  of  C,  may  be  constructed  by 
making  an  inner  projection  of  the  reaction  operator  t-  with  respect 
to  a  finite  basis,  where  the  reaction  operator  is  associated  with  the 
separation  of  the  Hamiltonian  Jf  ~  >f   +  V.   Several  types  of  inner 
projections  are  discussed  with  particular  emphasis  on  the  "LBwdin 
Projection." 

An  application  of  the  method  is  made  to  a  calculation  of 
lower  bounds  to  some  of  the  S  states  of  He  and  Li   using  the  "LBwdin 
Projection." 


CHAPTER  II 

THE  PARTITIONING  TECHNIQUE  AND  THE  BRACKETING  THEOREM 

Consider  the  problem  of  solving  the  SchrBdinger  equation 
J>i|^  =  E\lf,  where  the  Hamiltonian  }\  is  assumed  to  be  hermitian,  i.e., 
0(  =  J^   and  bounded  below.   Let  (jD  =  [  ^^  J  be  a  complete  ortho- 
normal  set  of  functions.   As  an  immediate  consequence  there  is  a 
resolution  of  the  identity. 


I  =  E-l9><?^  . 


(2.1) 


where  the  sum  over  the  index  i  implies  summation  over  the  discrete 
indices  and  integration  over  the  continuous  ones. 

Let  us  define  a  self-adjoint  projection  operator  (J)  such  that 
it  projects  out  of  the  space  spanned  by  ^  a  subspace  (^^   of  order  g: 


u)    a   subspace    ^h 


0=  |ll9><9/  . 

From  this  definition  it  follows  that   ^  satisfies  the  relations 

Let  P  =  1  -  ^  be  a  projection  operator  associated  with  the  subspace 
(jO^^f    the  orthogonal  complement  to  the  subspace  (©  : 


(2.2) 


(2.3) 


f.=f-S 


Q.  • 


It  then  follows  that 


p*=  p   ^     ©''=  iS)    ,  OP^P©=^o. 


(2.4) 


(2.6) 


The  function  y     can  be  projected  onto  the  complete  set  CXr 
and  the  SchrBdinger  equation  can  be  rewritten  as 

Operating  from  the  left  first  by  C?,    then  by  P,  we  obtain  the  set  of 
equations 

(©i8©  +  (95^P)iP=   FOV^,  (2.7a) 

(PJ?0    +  PS«P)'!P=    £PV^.        (2-7^) 

From  (2.7b)  we  find  that 

^Lpce-k)p1p^-  ^'-'^ 

From  {2,k)   we  see  that,  for  an  arbitrary  number  a, 
Then  (2.8)  can  be  rewritten  in  the  more  general  form 

(PK©)©^=  t'<S>+  P(E-«)Pj  P"^,  (2-9) 


Solving   (2.9)    for   P\|/,   we  obtain  the   expression 


Pf=Z^-0iPCE-5?)PTlPifS>)(^V.       (^-1°) 


Operating  on  (2.10)  from  the  left  by  P  gives 


(2.12) 


Substituting  (2.11)  for  P\|;  in  equation  (2.7a),  we  obtain 

(Oii<S>)&Pi-CS>iliP)P'iP  =^   BS>P 

(2.1 

= f  (9sf  (s>  f  (fi»i?p;  pr«<<5 + p(£- »)pr  w©]  ©v^ 

Let  us  introduce  the  reduced  resolvent  T„ 

-]j  =  prw.(p+p(£->f)pj'p .        (2-13) 

Equation  (2.12)  then  becomes 

Introducing  (2.11)  and  (2.I3).  into  (2.5)  we  obtain  the  expression 

P 
Let  us  introduce  the  symbolic  notation  T„  =  = rr  .   It  is  easily 

seen  that  T  satisfies  the  relations  P(E  -i^)Tp  =  P, 


Since  Orp^    ^'9!K5>JV^  -   ^  ^i  J^i    #     multiplying 
(2.14)  from  the  left  by  each  of  the  m   ^€i   C^^,  (k=l,2,  . . .  ,g)  and 
integrating  we  obtain  a  set  of  g  simultaneous  equations 


(k=l,2,...,g). 


(2.17) 


This  set  of  equations  can  be  rewritten  as  a  matrix  equation 

and  {C      is  a  column  vector  composed  of  the  coefficients  appearing  in 
the  expansion  of  C'r  in  terms  of  the  functions  ©.€•  ^P  . 

The  solutions  to  this  set  of  equations  are  found  by  setting 

the  determinant  of  the  coefficients  of  (L   in  the  matrix  equation 
—  »  a 

C  W^fc^^^  *  EJ^ai  ^a  ~  ^     equal  to  zero: 

The  solutions  of  the  secular  equation  (2.20)  are  eigenvalues  of  ^ 
associated  with  states  not  orthogonal  to  the  subspace  dp  .   This  point 
will  be  discussed  in  greater  detail  later. 

One  way  of  solving  equation  (2.20)  is  by  means  of  some  iterative 
procedure.   Let  £  and   c,  be  real  variables.  By  replacing  the 
variable  E  in  QrO       with  f  ,  the  elements  of  the  matrix  f/y       become 
functions  of  ^  and  the  secular  equation  (2.20)  has  a  set  of  g  solu- 
tions I  C 1 • f  ^°^  each  value  of  ^  ,  corresponding  to  those  values  of 


the  variable  £,  that  satisfy  the  equation 

Hence,  the  variable  £  ^  is  a  multivalued  function  of  £ 
having  the  same  multiplicity  as  the  order  of  the  subspace  (TD   and 
this  is  denoted  by  the  expression  C  ,  =  jT  ( £ ) • 

Corresponding  to  the  set  of  g  solutions  j  wi-j  obtained  by 
solving  (2.21)  for  a  fixed  value  of  ^  ,  there  will  be  a  set  of  g 
vectors  JsZ  a.^  . 

Making  the  substitution  of  £  for  E  in  (2.I9)  we  obtain  the 
expression 


where 


(2.23) 


Then,  for  arbitrary  £  ,  each  of  the  solutions  ^, .  ( i=l,2, . . . ,g)  may 
be  expressed  as  a  function  of£,  i.e.,  £^.   =f.(C),  having  the  form: 
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Hi  -     V  £  . —       =  Ti  ^^)  •  ^'-^'^ 

Taking  the  first  derivative  of  (2.2U)  with  respect  to  £.   and  evaluating 
it  at  the  same  value  of  £  for  which  the  (S  .  were  evaluated,  we  have' 


ai 


fiTlk'^^     A„    rf.  .  (2.25) 


We  see  from  (2.23)  that 


^'=  <gpJ.4(A)lf.)=-<fJTeV.v<-) 


We  then  see  from  (2.l6)  and  (2.23)  that 


from  which  it  follows  that 


Replacing  — — in  (2.25)  by  the  right  hand  side  of  (2.27). 

d£ 


we  obtain  the  result 


<  0     ,  (i=l,2,...,g) 

Therefore,  each  of  the  solutions  £ , .  is  a  function  of  g  having  a 
negative  slope  for  all  £  . 

Let  ^  =  E  +  6  ,  where  e  is  the  difference  between  the  value 
of  £  and  some  eigenvalue  E  of  V^  .   Let  one  of  the  solutions  given 
by  (2.2U)  be  represented  as   £    =  E  +  6,..   If  fe  is  a  small  quan- 
tity, we  can  use  the  Lagrange  mean-value  theorem  to  obtain  the  result 
that 


=  -^-(.B)  ■)-  &  ^'te-tfie)  ^ 


where  0  <,  0  ^    1.   Using  the  identity  E  =  f.(E)  (see  (2.18))  we 
obtain  the  relation 

From  (2.28)  we  know  that  the  first  derivative  f .  is  negative;  hence, 
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the  quantities  £  and  S^.    have  different  signs.   This  leads  to  the 

bracketing  theorem   which  says  that  the  pair  (f  »£•,.)  bracket  at 

least  one  true  eigenvalue  of  J^  .   Then  the  g  intervals  (£,  £,-,.)   will 

each  contain  at  least  one  true  eigenvalue  E  such  that  if  £  •  <  E, 

then  E  1 ■  y     E  and  vice  versa. 

**  1l  '  ____^ 

Let  us  examine  the  behavior  of  (rO     (£)  in  the  limit  of 

3d 

£  -*  +  GO.   We  see  from  (2.22)  that 

=  ^.»  +  M«.(  £^,.  <3>J  i?«  If,)}  ^6a   (-30) 

=  W)*a  • 

Hence,  the  set  of  g  solutions  {  fi-j   have  as  horizontal  asymptotes 
the  eigenvalues  of  the  matrix  0^      : 

Let  us  define  the  "outer  projection"  j/f    by  the  relation 

k=  PHP. 

Since  P©  =  ©P  =  0,  it  follows  that  all  the  functions  (D  .  €   (D . 
are  eigenfunctions  of  S^  associated  with  the  eigenvalue  0. 


(2.31) 
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Let  us  consider  the  normalized  eigenfunctions   ^.    of  [^ 
associated  with  the  eigenvalues  E.  that  are  situated  in  the  subspace 
(S)  .  We  then  have  the  spectral  resolutions 

p=  Eivjx^i  . 

*-  (2.32) 

It  is  obvious  that  the  reaction  operator  has  vertical  asymptotes  for 
the  values  ^  =  E.,  (i=l,2, ...)•   It  then  follows  that  the  bracketing 
function^  ^  has  vertical  asymptotes  for  the  values  £,  =  E.,  provided 
E.  is  not  simultaneously  an  eigenvalue  of  J*\  . 

Let   i|/.  be  a  function  satisfying  the  SchrBdinger  equation 
\Pilr.  =  E.ijf..  From  (2.5),  we  see  that  i|f.  =  (P^.  +  P\|f.  .  Therefore,  if 
(J)\lf.  =  0  and  P\|;.  ?^  0,  it  follows  that 

}p.     ^    P  lf>.    =     Tp.        ,  (2.33) 

Then,  from  (2.7a)  and  (2.7b)  we  find  that 


(2-31t) 


^^.  =   (P«P)^i  =  Eji/^. 

That  is,  for  those  values  £  =  E.   where  E.  is  an  eigenvalue  of  J^ 
associated  with  a  state  \|^.  that  is  orthogonal  to  the  reference  mani- 
fold O  ,   i|a.  is  also  an  eigenstate  of  J#  having  the  same  eigenvalue, 

i.e.,   E.  =  E. . 
1    1 

Let  us  make  the  separation 

P=  P'+  P".  (2.35) 
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where  P'  =  Ill^^K'^fJ      for  all  the  f^   =  t^^  such  that  (tf^if^  =  0  a 
a 

P"  =     2lltB)<^^J  ^°^   3^^  the  t   such  that  ©iL  j'  0. 

The  reaction  operator  Tg  may  then  be  separated: 


iS    C*    :i^         (2'.36) 


=         'e  +         U    . 

It  then  follows  from  (2.22)  that 

.  u    .  ri$LaiMiMiM_      (3.37) 

Since  j^  ^Qj  =  E^^c^  and  ©  t^^  =  0,  we  see  that  i(f^\)K\  ^^     = 

E  <©  I  t  /^  =  0  foi^  all  (p     e.  Cp .      Hence  the  second  term  in  (2.3?) 
drops  out  and  we  obtain  the  final  result 

where 

•From  (2.38)  it  is  easily  seen  that  for  S    =  E^,  where  E^  is  a 
solution  of  If  t.  =  E\|f.  and  G)t-  =  0»  the  bracketing  functionf^^  has  no 
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vertical  asymptote.   Further,  we  see  that  S  ■,    has  vertical  asymptotes 
only  for  those  values  of  £  corresponding  to  eigenvalues  of  P  ^  P  not 
simultaneously  eigenvalues  of  J^ ,  i.e.,  only  for  those  cases  where  (P 

^i  ^   °' 

o 

It  follows  from  the  variation  principle   that  the  eigenvalues 
of  ^  are  upper  bounds  to  the  eigenvalues  of  ff   in  order, 

£^  >/   E  J  .  (2.39) 

Hence  the  vertical  asymptotes  appear  in  the  bracketing  function  £ 

for  values  of  £  that  are  upper  bounds  to  eigenvalues  of  /f  corresponding 

to  the  states  \|/.  satisfying  the  criterion  that  ^  \|f.  '^   0. 

A  plot  of  the  g-valued  function  £  as  a  function  of  the  vari- 
able £  can  be  made  by  imposing  the  non-crossing  rule,  which  is  valid 
for  states  having  the  same  symmetry.  The  values  of  the  g  curves  will 
satisfy  the  relation 

£fl     >/      Ef  ,     (i=l,2,...,g),      (2.i+0) 

for  all  £  ^  E^ .   If   £^^  is  taken  to  be  the  lowest  value  of  the  set 
of  solutions  j^-,.T    ,  then  for  £  chosen  equal  to  E^  we  find  that 

C--,  =  E   provided  that  <D\|/,  ^  0.   Then  for   £  >  E^^,  we  see  from 
(2.29)  that   c  IT  ^  E^  and  the  C-,-i    curve  will  undergo  a  discontin- 
uous change  from  -  00  to  +  00  as  £  passes  through  E^ .   In  general, 
we  find  that  the  g  curves  have  the  behavior  illustrated  in  Fig.  1, 
which  is  for  g  =  3*   We  see  that  each  branch  of  the  curves  crosses  the 

ti  =  z.    line  at  an  eigenvalue  of  ^  associated  with  a  state  \|f.  not 
orthogonal  to  the  reference  manifold  (©  . 
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Fig.  1.  Typical  behavior  of  the  multivalued  bracketing 
function  £'j^  =  f (  E  )  f or  g  =  3. 
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It  is  apparent  from  this  discussion  that  for  ^  ^.    ^  0   and 

^-  K  S  K    ^•»  O"^^  of  ^^^   eigenvalues  of  (}v  (6)>  C  ii,>  i-s  the  lower 
bound  to  E.,  i.e.,  E.  lies  in  the  interval   (f  •,,,£)  and  for  f  =  E., 

*"  Ik    1 


CHAPTER  III 


LOWER  BOUNDS  TO  ENERGY  STATES 


Assume  that  the  Hamiltonian  1/  of  a  system  may  be  separated 
into  two  parts;  a  zeroth  order  Hamiltonian  or  unperturbed  part  ^      and 
the  interaction  part  or  perturbation  V.  Let  there  be  associated  with 
the  Hamiltonian  i4      a  complete  orthonormal  set  of  eigenf unctions  i  Q) .( 
and  a  known  eigenvalue  spectrum  E.,  (i=l,2, ...)»  such  that 

Let  it  be  further  assumed  that  the  perturbation  V  is  positive  definite: 

V  >  0   .  (3-2) 

The  space  (p         spanned  by  the  set  1 CP  • (  can  be  partitioned 
into  two  subspaces  (|p    and  (^  ,  by  defining  the  projection  operators 

©=ii0<<3^i  =  t-o: ,  (3.3) 


L«/  ^  t»/ 


p=  1-0=  ii(?f><9;'  -  ^^i ' 

such  that  '   <'  V    d 


The  operators  0  and  P  satisfy  the  relations  given  in  (2.3)  and  (2.U) 
From  (3'l)  we  also  see  that 
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and  using  (2.l6)  we  have  the  result 
Equation  (2.I9)  may  be  rewritten  as 

Let  us  introduce  a  reduced  resolvent  T  associated  with  the  unperturbed 

o 

Hamiltonian  X  , 

where  again  the  sum  implies  summation  over  the  discrete  states  and 
integration  over  the  continuous  part  of  the  spectrum. 

For  any  inverse  operator  of  the  type  (A-B)   there  are  the 
following  identities: 

Lfl-Sr'  =   B''  +  fl''BiB-B)'! 
=  /?■'  4-(^-3)"'6^    . 

provided  the  inverse  operators  involved  actually  exist.   Using  the  first 
of  the  identities  given  in  (3.7)  we  find  that 

*  (3.8) 
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where      :t'^   ^     V+VT^V     £       V     +     V  X"  jfg  (3-9) 

2 
is  the  reaction  operator.   The  last  relation  corresponds  to  the 

Lippmann-Schwinger  integral  equation  in  scattering  theory. 

Provided  V   exists,  we  find  from  (3.9)  that  (V~  -T  )tg  =  1; 

hence, 

A  =  iv-xy'-  <5-^°' 

Assuming  the  spectrum  of  ^f   is  ordered,  E^  <  E  <  E  <  ..., 
it  is  easily  seen  from  (3.6)  that  for  f  <E  ^,  T   <  0.   Since  it 
is  assumed  that  V  >  0,  we  see  that  V   >  0  and  hence  t^  >  0. 

Introducing  (3'9)  into  equation  (3«5)  we  obtain 

Hence,  we  see  that  the  choice  of  ^  as  given  in  (3-3)  leads  to  a 
multi-valued  bracketing  function  ^  ^ .   From  Chapter  II,  it  foil 
that  this  bracketing  function  will  have  vertical  asymptotes  only  for 
those  states  \|f,  of  jt    that  satisfy  the  relation  C)i|^^^  f'  0.   That  is, 
only  those  states  not  orthogonal  to  the  reference  manifold  (ff)   , 
spanned  by  the  set  of  g  eigenvectors  of  V  ,  will  cause  vertical  • 
asymptotes  to  appear  in  ^  ■,  -      The  function  ^^    has  tne  important 
property  that  the  first  vertical  asymptote  occurs  for  £  >  E  , ^ . 
Hence,  if  P  Jf  ^  has  the  ordered  spectrum  E^  <  E  <  ...,  then  E^  >  E   ^. 
To  prove  this  let  us  consider  two  self -adjoint  operators  A  and  B 
satisfying  the  relation 

/)      >       B.  (3.12) 

Then,  for  any  function  y,   <  O  / /)  I  ^>   >  ^<flBl<^'>' 


(3.11) 


Lows 
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Let  P  be  an  arbitrary  linear  operator.   Then  replacing  (p   by  P^  we 
find  that  <q)|P+APl  (5!>>^(plP+BP|»>   ,  i.e., 

P'^RP  >    P^BP.  ^'-''^ 

Since  V  ;>  0  we  have  the  inequality  i^  =  if  "^  ^  ^  H  '  ^'^^^^  iC*"  ~  W 
ar.d  ^/>  ~  =  J?  .  Then  from  (3. 13)  we  find  that  for  the  self-adjoint 
projection  operator  P  defined  in  (3-3) 

Hence,  <0  |P^PJ0)><^^l  Pi{<^|(p^    for  any  function  ep.   Let  ^  =  ijr^, 
where  y,  is  the  lowest  state  of  P|;fP  not  orthogonal  to  (^     .   Since  ■ 
the  ground  state  of  P^^P  is  E  .  ■,  >  it  follows  from  the  variation  princi- 
pie   that 


It  follows  at  once  that  the  first  vertical  asymptote  of  cu   appears 

-54        O 

for  '^^  >  E  ,  1  •   In  general,  the  eigenvalues  of  a  self -adjoint  operator  A 

are  greater  than  or  equal  to  the  eigenvalues  of  a  self-adjoint  operator 

2 


.  B  in_  crder  if  A  >  3 . 


It  is  an  immediate  consequence  of  this  that  for  £  <    E   ,  tg_ 


2,9 
Lsts  and  s^^isfies  the  inequality  ^ 


Henc 


find 


tg^    -      ^''-"77      ^    ^"'     >     ^-  (3.16) 

e  V  ?  tj    >   0  and        Kq  "^  ^  ^     ^o  "^  ''^  ^  ^o'      ^PP^^^^S   (3-13)  we 
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and  therefore. 


-^.      >    £f-     >  Bl^      (i=i,2.....g).      (3.17) 

This  says  that  the  eigenvalues  of  ^  ^  are  bounded  above  by  the  average 
values  of  ^   with  respect  to  the  basis  (SD      and  bounded  below  by  the 


eigenvalues  of  vr  ^il  order  for  £  ^  E  , ,  . 


g+1 
For  each  value  of  c  »  C,  will  satisfy  the  equation 

Lt,-E,j)   £;   =0,         ('-^'^ ^)-    (3.18) 

From  this  it  follows  that  for  any  ^  the  curves   Ci-»  (i'^ljS,  . .  •  ,g)> 
have  negative  definite  slopes.  From  (3.IO)  and  (3.II)  we  have 

(3.19) 

dT       2 
VJe  see  from  (2.26)  that  -  — «-  =  T   .   Making  this  substitution  in 

(3.19)  we  obtain 

Each  of  the  curves  has  a  slope  at  the  point  £  for  which  ^  ,^   and  (J 
were  evaluated  given  by 

dZ,[  I        ^    ct  «J?'  (t-  f        <  0  ,      N 

(  ial>2, .  . . ,g)  . 
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Hence,  the  curves  have  negative  definite  slopes  and  the  bracketing 
theorem  holds  for  each  of  the  g  eigenvalues  of  £  . .   From  (3.2I)  and 


the  inequality  (3.I5)  it  follows  that  for  every  E^^  <  ^  <  E^j^»  where 
E,  is  the  k  th  lowest  eigenvalue  of  the  Hamiltonian  ^  associated  with 
a  state  \|;,  not  orthogonal  to  the  subspace  ^)  ,  there  is  at  least  one 
solution  to  (3.18)  that  will  provide  a  lower  bound  to  E,  .   In  fact,  it 
is  now  easy  to  prove  that  under  these  conditions  the  k  th  lowest 
solution  of  (3.18)  is  the  lower  bound  (or  highest  lower  bound)  to  E,  . 

From  the  discussion  associated  with  Fig.  1,  we  observe  that  the 
curves  of  ^  ■,    cross  the  Ci  =S  Hfie  at  the  eigenvalues  of  >^  asso- 
ciated with  states  f.    having  the  property  that  (5)\tf.  ^  0.     From  (3*  15) 
we  found  that  the  first  vertical  asymptote  was  situated  to  the  right  of 

E  , T ,  i.e.,  Et  >  E  , T .   Let  there  be  k  eigenvalues  of  If  having  the 
g+i         i    g+i  '» 

ordering  E^  <  E  <  ...  <  E  <  E  ^,  associated  with  states  satisfying 
the  above  property.   Let  the  eigenvalues  of  ^  -,    also  be  ordered,  i.e., 
£^,,  -  £',p  <  •••  <f  1  •   Then,   f  ^^  is  the  lower  bound  to  E^  for 
E^  <  S  <C  E  ,  ^  and  is  equal  to  E^  for  f  =  E  .   In  order  to  see  this 
.let  us  assume  there  exists  for  some  j  >  1  a  solution  ^  .  such  that 
C 1 ^  ~   E,  for   C  ~  E^.   Then,  from  the  non-crossing  rule  and  (3.2I), 
£ -.-^   must  be  <  E,  and  hence  crosses  the  S ■,    =   C  line  for  some  value  of 
iJ  <  E^.   This  contradicts  the  assumption  that  E^  is  the  lowest  eigen- 
value of  y^    associated  with  a  state  f   satisfying  the  requirement  that 
Q^   ^  0.   Similarly,  Cip  provides  the  lower  bound  (or  highest  lower 
bound)  to  E  for  E  <  £  <E    .      If  for  some  j>  2   there  exists  a 
solution  £       such  that  S  i-   ~  ^p  ^°^      2  ~   ^  ,  then  by  once  again 
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imposing  the  non-crossing  rule  and  (3.2I)  we  find  that  tip  ^^^st  cross 
the  £*  ,  =2   line  for  a  value  E  <  ^  <  E  which  leads  to  a  contra- 
diction.  By  induction,  we  see  that,  for  any  E.  <  £  <  E  ^^,  Sn   ^^ 
the  lower  bound  to  E.,  ( i=l,2, . • . ,k) . 

Figure  2  illustrates  the  typical  behavior  of  the  multivalued 
bracketing  function  (3.II)  for  g  =  i^  and  k  =  2. 
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Fig.  2.   Typical  behavior  of  the  multivalued  bracketing 
function  \p      for  g  =  U  and  k  =  2. 


CHAPTER  IV 


LOWER  BOUNDS  TO  THE  REACTION  OPERATOR 


In  principle,  equation  (3.II)  can  be  solved  exactly  and  provides 
lower  bounds  to  those  states  of  the  Hamiltonian  iP    satisfying  the 
criteria  that  ^i|/^  f   0  and  £  <  E°   .   In  fact,  there  are  few  prob- 
lems for  which  (3.II)  can  be  solved  exactly. 

From  (3.10)  we  see  that  the  reduced  resolvent  T  must  first  be 

o 

obtained  in  order  to  calculate  the  reaction  operator  tg  .   From  (3.6) 

we  know  that  T  <  0  for  £   <  E°  . .   If  T  cannot  be  evaluated  exactly, 
o  g+1       o 

one  approach  would  be  to  estimate  it  by  the  operator  T  (p): 


o' 


E. 

(it.i) 


where  it  is  assumed  that  the  eigenvalues  are  arranged  in  the  order 
E,  <  E  <  ....   T  (p)  is  an  upper  bound  to  T  .   Hence,  from  (3.IO)  we 


see  that 


and  since  it  is  assumed  that  V  >  0, 


2k 
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-/       .       /  -J 


Since  tp  (py  and  tg  are  positive  definite,  the  inverse  operators  t£^(p) 


2.9 
and  tc        exist,  are  positive  definite   ^  and  satisfy  the  inequality 


t^    >y    :&^Cf)     >  0  . 


{h.k) 


Better  lower  bounds  to  the  reaction  operator  maybe  gotten  by 

introducing  the  idea  of  inner  projections. 

Let  f  =  (fT,f^,...,f  )  be  a  set  of  n  linearly  independent 
^  1  2      n'  J  r 

vectors  spanning  a  subspace  in  Hilbert  space  having  the  metric  A  =  <f J t^. 
Then  the  projection  operator  associated  with  this  linear  manifold  is 
given  by 

4J 


Q=il>A-'<|)--  i.ii>h:Uij.   (...) 


Let  P  =  1  -  Q  be  the  projection  operator  associated  with  the  orthogonal 
complement  to  f. 

The  projection  operator  Q  is  positive  definite  since 

Hence,  P  =  1  -  Q  ^  0  since  P  also  is  a  projection  operator.  Therefore, 
Q  satisfies  the  inequalities 


b^'  q>^i 


(h.i) 
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The  reaction  operator  t^  is  positive  definite  for  P  <E  ^,  and  self- 

«-  *•    g+l 

+  -  11 

adjoint:  ^i   ~  ^£,   •      Hence,  t^  exists  and  is  positive  definite.    The 


inner  projection  of  t^  with  respect  to  the  subspace  f  is  given  by 


(h.Q) 


(i^.io) 


From  (U.7)  we  find  that  tl  satisfies  the  inequalities 
Let  us  define  two  new  manifolds  ^  and  ^  ,  where 

Substitution  of  these  expressions  into  (4.8)  results  in  the  following 
relations : 

and 

The  inner  projection  given  by  (U.ll)  will  be  called  the  Aronszajn  projec- 
tion and  the  space  (9  =  (g, ,gp, ...,g  ),  the  Aronszajn  space.   That  given 
by  (^4^.12)  will  be  called  the  Bazley  projection  and  the  associated  space 
^(k    =  (h^,h  ,  ...,h  ),  the  Bazley  space. 

LBwdin  has  shown  that  a  lower  bound  to  an  energy  level  E  of  the 
Hamiltonian  W  may  be  obtained  when  there  are  g  eigenvalues  of  vP 
less  than  E,  if  the  Aronszajn  space  (&  contains  at  least  all  the 
functions  ^i»(&p»  •••»©   associated  with  the  g  unperturbed  levels 
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and  the  reference  function  <D   is  constructed  by  taking  a  linear  combi- 
nation of  these  g  functions. 

Replacing  the  reaction  operator  t^  in  (3-ll)  with  the  Bazley 
projection  defined  by  (ll-.12),  we  have  a  new  bracketing  function 


£^=<fMjttif:y. 


(iv.13) 


From  {k.9)   we  see  that  for  S   ^   ^g+i' 


^At,  >.  ^^+t^  >  Z  ("-^^^ 


and  therefore 


From  this  we  may  conclude  that  the  eigenvalues  of  ^,  are  bounded  above 
by  the  eigenvalues  of  ^,  and  bounded  below  by  the  eigenvalues  of  W 
in  order; 

From  the  bracketing  theorem,  g^.  is  a  lower  bound  to  E.,  where 
E.  is  the  i  th  lowest  eigenvalue  of  |f  satisfying  the  criteria  E  ,  j^>  £  >E. 
and  £P  y.  f'  0.   Therefore,  ^'   is  also  a  lower  bound  to  the  same  eigen- 
value.  From  (3.17)  we  observe  that  the  eigenvalues  of  *|  ^  are  bounded 
above  by  the  eigenvalues  of  ^  with  respect  to  the  basis  ®   in  order. 
It  follows  at  once  that  the  latter  also  provide  upper  limits  to  the 
eigenvalues  of  ^,  in  order. 


28 


More  favorable  upper  bounds  to  the  eigenvalues  of  ^^  may  be 
obtained  by  introducing  the  inequality  (3.I6)  into  the  expression  for 
tg  which  gives  the  result 

Hence,  we  have  that 

where  ^im  t^   =  V  and  V  is  simply  a  Bazley  projection  of  V  onto 
the  space  spanned  by  the  manifold  vV.  It  follows  at  once  that 

V  >   V^  >  0  .  /  (U.18) 

From  {k.lT)   we  see  that 

hence. 

The  eigenvalues  of  ^  ^    are  bounded  above  by  the  eigenvalues  of  ff   +  V 
and  bounded  below  by  the  eigenvalues  of  Jf  : 

hi     ^   £,i    9-  Bl    ,  (i=i,2,....g),  (iv.2o) 

where  h.  is  the  i  th  lowest  eigenvalue  of  V\  with  respect  to  the 
basis  (S>  .   The  inequalities  given  in  (i4-.20)  hold  for  all  ^  <  E    and 
tell  us  at  once  that  ^ '  has  no  vertical  asymptotes  in  this  region  for 
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any  Bazley  projection.  The  equality  £^.   =   E.  occurs  whenever^^  =  0: 
when  the  Bazley  space  is  chosen  orthogonal  to  the  reference  manifold 

The  inner  projected  bracketing  function  given  in  (U.I3)  can  be 
computed  exactly  in  only  a  few  cases  involving  particular  choices  of  the 
Bazley  space;   '  ~^'    otherwise,  T  is  truncated  according  to  (li-.l), 
which  results  in  an  operator  tg  (p)  that,  for   ^  <  E   ^,  satisfies  the 
inequalities 


^E  ^  ^s"  >:^g^c^>    >o    . 


(1^.21) 


To  remove  the  problem  of  dealing  with  the  reduced  resolvent  T  , 

o 

Lbwdin  has  suggested  the  substitution  /j  =  ( C -  Jf  )>f  be  used  in 

making  the  inner  projection.   Substituting  this  into  (14-.  12)  we  find  that 

(1^.22) 


where 


since,  from  (2.I6).  (  £  -  Sf„)T^  =  (  £  -  ffo)PT^  =  P(  £  "  i{^)'r^  =   P. 
We  shall  call  tg  defined  by  (U.22)  the  LBwdin  projection  and  the  asso- 
ciated space  A    =  ( jj^,  jp,  .  . .,  j^),  the  Lbwdin  space.   This  projection 

has  the  interesting  property  that  the  reduced  resolvent  T  has  been 

o 
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eliminated  from  the  calculation.   It  also  leads  to  some  very  interesting 
behavior  not  observed  wnen  using  the  Bazley  or  Aronszajn  projections. 
This  behavior  can  be  explained  by  the  fact  that  the  linear  manifold  onto 
which  the  projection  of  tg.      is  made  is  itself  a  function  of  the  variable 
^  and  varies  with  changes  in  ^ . 

An  immediate  consequence  of  the  C  dependence  of  the  Lbwdin 
space  is  that  the  slopes  of  the  eigenvalue  curves  of  ^,  are  no  longer 
negative  definite.   In  fact,  minima  occur   in  these  curves  for  values 
of  t  =  E.,  (i=l,2,  .  . .  ,g) .   That  is,  one  of  the  eigenvalues  of**,  for 
C  -   E,  ,  E,  <  E  ,  ^  and  ^i,^<|P  >  ist,^  =  E,  and  E,  will  be  a  minimum 
value  of  the  function  C  io* 

This  is  easily  seen  when  we  write  down  the  expression  for  the 
iJl   th  element  of^,  using  the  LBwdin  projection  (i4-.22). 


The  only  non-vanishing  term  in  the  k  th  row  or  k  th  column  of ^ ,  when 
I;  =  E°  is  the  kk  element,   (Cl^kW  ~  ^k'  ^^   ^^'^   write  the  expression 
for  (U.I5)  evaluated  at  the  point  C  ~  E,  , 
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Vii  ••• 


•lk-1 


t' 


•  •  •      w« 


i[ 


i-< 


■k-1,1 
0 

.1 

'k+1, 1 


^_,+t,'    ,    „  ,      0 


■1     k-l,k-l 
0 


-k+l,k-l 


*^k-l,k+l         •••   ^-1,1 


Q 


0       Ef,,+t; 


k+1   "k+l,fcfl 

I 


Q 


'k+l.g 


'•8,1 


-g,k-l 


'g,k+l 


E°+t' 
g     gg, 


(,k.2k) 


Obviously,  E,  is  an  eigenvalue  of  ^'    for   ^=   E".   That  this  is  also  a 


^k- 


minimum  in  some  eigenvalue  curve,  e.g.,  p\p  °^  £ -i  Is  seen  by  looking 
at  the  first  and  second  derivatives  of  ^  -,  with  respect  to  ^  evaluated 
at  C=   E°.   If  £1^  =  E°  =£  ,  then 


-  o 


£--^ 


>  o 


are  sufficient  criteria  for  the  point   £  =     E   to  be  a  minimum  in  the 
&  In   curve.   From  (3.21)  we  know  that 


d£.:  =  c/  e^' 


de 


4£ 


&. 


6.e| 


Ct.a?) 
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where  (J  »  is  thei  normalized  eigenvector  of  (T  ^  associated  with  the 


eigenvalue  E,  .  Since 


k\ 


0 


is  the  eigenvector  for  £  ip  "^  '^i-  '^  ^ '    ^^   follows  from  (14-.25)  that 


(0  ...  01,  0  ...  o) 


£*Et 


£.6° 


dt, 


dz 


a 


£•£■ 


(It.  26) 


.-/ 


From  (U.23)  and  the  identity  (i  A   S  - 


4"'il  A'' 


.  we  obtain  the  expression 


\d€ 


=  [2i€-^<<p;i^>A-'<0:>- 


(U.27) 


from  which  it  follows  that 


c?2 


6^e; 


=     6>. 


(1^.28) 


£^Bl 
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The  expression  for  the  second  derivative  is  given  by 


ST       /rt     Of       S, 


8=E 


de"- 


£*E 


■^ 


=  [^  <<pri^>  ^'"<^l  9P  ~  "^  (^-^P  <5^li^> 


OvA 


'■*"'^<^¥>-<e-2.")'<9.*'^> 


^{A-'#/S-'i<^i?J)] 


£r=e"^ 


;i<(p;/^>A"'<^lfC>     >o. 


(i^.29) 


Hence,  from  (l|-.28)  and  (k.2^)    it  is  obvious  that  the  g  curves  of^^ 
using  the  Lliwdin  projection  are  no  longer  monotonically  descending. 

In  general,  the  curve  £|  •,  «  can  have  minima  and  cross  the  £  !.    =£ 
line  at  several  of  the  points  £)  =  E..   It  may  also  cross  the  £.    =£ 
line  at  values  of  ^    not  eigenvalues  of  ^  .   From  (l<-.20)  we  have  the 
relation   h.  ^  £    .  >y    E.    for  any  £  <  E°   .   It  follows  that  £,     . 
can  have  minima  at  any  or  all  the  points   £  =  E.,  where  i  ^  X  and 
E.  4h^.   These  features  of  the  LBwdin  projection  are  developed  further 
in  the  application  of  thi-^  theory  to  He  and  Li  . 

In  this  section,  we  have  developed  a  procedure  for  calculating 
lower  bounds  to  the  eigenvalues  of  a  Hamiltonian  t\  =  CP  +  V,  where 


3i^ 


V  >  0,  and   \y      has  a  set  of  known  eigenfunctions .   The  subspace  (^ 

is  chosen  such  that  ^  if .    "^  0,   where  \j/.  is  the  eigenf  unction  of  ^ 

associated  with  the  i  th  lox^/est  eigenvalue  E.  having  a  particular 

symmetry  and  such  that  it  contains  all  the  eigenfunctions  associated 

with  eigenvalues  of   vP   less  than  E.  having  the  same  symmetry  as  f'.. 

The  L'dwdin  space  -i  must  be  chosen  such  that  the  overlap  integrals 
or  F     & 


\(p-|j,^  with  the  functions  of  the  reference  manifold^?  are 


re  non- 
zero.  It  is  then  a  rather  straightforward  matter  to  calculate  the 


matrix  elements  given  in  (i|-,23).   If  1'i^((^)  =  gj  then  by  choosing  the 

parameter  p  >  E.  and   <E  ,^,    the  i  th  lowest  solution  of  the  secular 
«-    1         g+1 

equation 

orovides  a  lower  bound  to  E.. 

X 


CHAPTER  V 


LOWER  BOUNDS  TO  THE  "'S  STATES  OF  He  AND  Li 


The  procedure  discussed  in  Chapter  IV  is  used  to  calculate 

o  + 

lower  bounds  to  some  of  the  low-lying  S  states  of  He  and  Li  .   For  the 

two-electron  series,  the  nonrelativistic  Hamiltonian  in  atomic  units 
is  given  by  "      • 

/ 


(5.1) 


where  z  is  the  nuclear  charge,  r,  and  r  are  the  position  vectors  of 
the  electrons,  r .  =  |r . | ,  i  =  1,  2  and  r^  =  j r^-r  j . 

We  will  take  advantage  of  the  symmetry  properties  of  the 

o 
Hamiltonian  (5.I)  and  restrict  ourselves  to  states  having  S  sjanmetry. 

For  these  states,  the  eigenf unctions  of  (5'l)  depend  only  on  r,,  r  and 

r,   and  are  antisymmetric  in  the  spatial  coordinates  of  the  two  elec- 


trons, that   is 


RzWcr. ,^)  =  - '^Cf. ,n)    , 


where  P^   is  the  permutation  operator  acting  on  the  coordinates  r^  and  r  . 
The  Hamiltonian  (5.I)  may  be  separated  into  an  unperturbed  part 
^*r   and  a  positive  definite  perturbation  V,  where 

i^o"^   '''^'^^"  ^"^^  ""it  ""i^        '        (5.3) 

1/  =  -fc  >  0  . 
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The  eigenf unctions  of  I,"   are  well  known  and  for  the  S 

6\   o 

symmetry  are  given  by 


where  the  R  „   are  normalized  hydrogen  radial  wave  functions,  P,,  is 

the  £    th  normalized  Lengendre  polynomial  and  ©^   is  the  angle  between 

the  radial  vectors  r^  and  r^.   The  C    „  are  chosen  in  such  a  way  that 

1      2        '^I'^p-^ 
o  2 

the /K>     /,  are  normalized  to  it  , 

-»   1^2 

,/     /©\>    =  TT^.  (5.5) 

The  eigenvalue  spectrum  for  the   S  levels  of  J^_^  is  given  by 


I:  [  <  "^  <l  ■  '^'■^^ 


For  the  lowest  ordered  S  states  of   M*   ,  one  electron  is  in  the  n  =  1 

2 
level  and  the  angular  momentum  L  has  the  eigenvalue  zero,  ^  =  0.   The 

eigenfunctions  and  eigenvalues  of  W      are  then "given  by 


9f  =5fi:  [f?.o<''>  l^nJ^'  -  Ri-^) /?..,^) J  , 


(5.7) 
)  . 


By  including  in  the  reference  manifold  {?y   only  those  functions  ^^. 
having  "^S  symmetry,  we  obtain  from  (ij-.30)  lower  bounds  only  to  those 
states  of  \v'  having  this  same  symmetry  since  ^^\]r   =  0  for  any  state  \|/ 
of  Cj.'not  having  S  symmetry.   Hence,  the  solutions  of  (^-30)  will  be  in 
order  lower  bounds  to  the  2'^S,  3  S,  l4--^S,  ...,  j  S  states  of  JT  ,    i.e.. 


37 


for£<E°^^  and£>E.3g, 

With  the  solutions  of  jr  available  for  the  construction  of  the 
reference  manifold^  ,  we  must  next  select  the  functions  for  the 
LBwdin  space  J^  .      It  has  already  been  pointed  out  in  Chapter  IV  that  the 
choice  of  functions  for  the  A   nianifold  should  be  made  subject  to  the 
condition  that  the  overlap  integrals  ^(D.jjj^  be  non-zero.   Funct: 
satisfying  the  criteria  that 


:ions 


P. 


'•xii^'i^  • 


(5.8) 


where  L  is  the  angular  momentum  operator,   can  satisfy  these  conditions. 

The  functions  were  selected  on  the  basis  of  those  introduced  by 

Is 
Hylleraas   and  are  given  by 

^.  =  c,^,.^^  e''*s*^  i-^'*'  u^'-  .  (5.9) 

where  y\     is  an  arbitrary  scaling  parameter  and 

S-Y,-<rY^      ,    t-  r,-T^     .    U-Y,^    .       (5.10) 

The  C,  „         are  normalization  constants  selected  in  such  a  way  that  the 

k.^.m. 
Ill 

o 
functions  given  by  (5.9)  ^re  normalized  to  l6it  .  This  choice  of  normal- 
ization proves  to  be  convenient  in  the  calculation.  It  is  easily  seen 
that  the  functions  given  by  (5«9)  satisfy  the  criteria  (5.8). 


38 


The  minimum  size  of  the  reference  manifold  is  determined  by  the 
number  of  states  of  iP  having  the  same  symmetry  as  and  eigenvalues 
less  than  those  of  the  state  of   Cr  to  which  a  lower  bound  is  sought. 

For  the  S  states  of  He  the  ordering  in  atomic  units  (a.u.)  is 

given  by 

E°  =  -2.500  ...  ^E°  =  -2.22  ...  <  E  _   =  -2.175229   '^^  < 
i  ^  2-^S 

E°  =  -2.1250  ...  <  E°  =  -2.080  ...  <E  ^  =  -2.0687  ^®  < 
3  ^      .        33s 

E°  =  -2.05050  ...  (5-11) 

and  for  Li  , 

E°  =  -5.6250  ...  <E  _  =  -5.110727  ^^  <E°  =  -5.00  ...  < 
1  23s  2 

E°  =  -4.78I250  ...  <E  .  =  -i+.752  ^^  <E°  =  -i+,68  ...  < 
3  33s  ^ 

E  3  =  -4^37  ^^  ■<■€   =  -4.6250  ...    .  (5.12) 


From  (5.11)  it  follows  that,  to  calculate  a  lower  bound  using  (4.30)  to 
the  2  S  state  of  He,  we  must  at  least  include  CP .    and  (P   in  the 
reference  manifold  (K'  »  i.e.,  g  =  2  and  for  a  lower  bound  to  the  3  S 
state  we  must  choose  g  =  k,    i.e.,  (S?      =      (D  i »  ^  o>  ^  o»  ^),*  Erom 


f:-  ^i- 3>^- 9;- 9 


h- 


(5.12)  we  see  that,  to  obtain  lower  bounds  to  the  2  S,  3  S  and  k  S   levels 
of  Li  ,  it  is  necessary  to  have  g  =  1,  3  ^^^   4  respectively.   Of  course, 
as  many  additional  functions  as  desired  may  be  included  but  these  values 
of  g  give  the  minimum  sizes  of  the  reference  manifolds  necessary  in 
order  to  calculate  lower  bounds  by  this  procedure  to  the  levels  indicated. 

The  size  of  the  Lbwdin  space,  which  shall  be  denoted  by  NINT, 
can  vary  from  one  to  all  the  functions  of  a  complete  set.  If  the  set 
is  compli^ce,  then  the  exact  eigenvalues  of  ^rx  are  found  to  be  those 
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values  of  q      at  which  the  curves  ^  ^  .    cross  the  c  i  ~  £  line.   In 
general,  however,  this  set  is  finite  and  the  curves  cross  the  C i  ~  C 
line  at  the  eigenvalues  of   VP   associated  with  functions  in  the 
reference  manifold.   Crossing  may  also  occur  for  values  of  fc  jf'  E., 

i  <g  +  1. 

Once  the  reference  and  internal  projection  manifolds  are  chosen, 
the  matrix  elements  of  F-,t    given  by  {k.23)t   Jnay  be  found  and  the 
secular  equation  (^.30)  solved.   Both  of  the  parameters  ^  and  ^  are 
varied  to  give  a  maximum  in  each  of  the  g  lower  bound  curves,  where  the 
variation  of  Sp!   with  respect  to  ^  is  subject  to  the  condition  that  £| 
be  greater  than  the  energy  of  the  state  to  which  a  lower  bound  is 
desired  and  less  than  E  . .   The  computation  procedures  involved  will  be 
discussed  in  the  Appendix. 

We  see  in  Table  1  the  lower  bounds  obtained  by  this  procedure 
using  1<0  functions  of  the  type  given  in  (5.8)  to  construct  the  LBwdin 
projection,  where  the  restrictions  k.  =  0,  1,  2,  Z.   =  0,  1,  2,  3,  k, 
m.  =0,  1,  2  and  k.  +  m.  ^3  were  imposed  on  the  powers  of  the  variables 

s,  t  and  u  respectively. 

■3 
Lower  bounds  were  calculated  for  the  2  S  level  of  He  by  taking 

the  size  of  the  reference  manifold  to  be  g  =  2  and  g  =  k,    and  for  the 

3  + 

2  S  level  of  Li  by  taking  g  =  1  and  g  =   k.      Since  the  size  of  the 

reference  manifold  determines  the  size  of  the  secular  equation,  a  fourth 

order  equation  is  the  largest  to  be  solved.   From  the  results  in 

Table  1  we  observe  that,  so  long  as  all  the  eigenfunctions  of  ^ 

associated  with  eigenvalues  less  than  the  2"^S  level  of  U>  are  included 
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in(^^°,  little  improvement  is  found  when  more  reference  functions  are 
added.   The  magnitudes  of  the  errors  in  the  bounds  to  the  higher  states 
are  seen  to  be  much  greater  than  that  of  the  2  S  state  in  both  the  He 
and  Li  calculations.  This  is,  however,  not  an  unexpected  result.   It 
has  been  pointed  out  by  LBwdin   that  a  truncated  set  which  works  well 
enough  for  the  ground  state  may  give  a  surprisingly  poor  result  for 

even  the  first  excited  state  of  the  same  symmetry. 

17 
Quite  good  lower  bounds  have  been  calculated  by  Pekeris   to  the 

o  20 

2  S  state  of  He  using  Temple's  method.    By  solving  a  determinant  of 

order  252,  he  obtained  a-  lower  bound  of  -2.175239  a.u.   The  lower  bound 

to  this  state,  for  g  '  k,   which  is  given  in  Table  1,  is  very  close  to 

this  value.   In  order  to  obtain  a  lower  bound  to  an  energy  level  of  ^ 

using  Temple's  formula,  it  is  necessary  to  calculate  matrix  elements  of 

2 
v#   and  to  have  a  knowledge  of  the  next  higher  level  of  the  system. 

It  is  apparent  from  the  work  of  Pekeris  that  in  order  to  obtain  a  good 

lower  bound  using  Temple's  formula,  the  size  of  the  secular  equation 

must  be  very  large. 

21 
Goscinski   has  shown  in  his  dissertation  that  Temple's  formula 

can  be  derived  from  (2.I9)  for  g  =  1. 

22  3 

Miller   has  obtained  lower  bounds  to  the  low-lying  S  states 

2' 
of  He  using  the  method  of  intermediate  Hamiltonians .  "'     Solving  a 

secular  equation  of  order  9>  he  calculated  lower  bounds  to  the  2  S  and 

3^S  levels  of  -2.l802  a.u.  and  -2.070i»-  a.u.  respectively.   His  result 

for  the  3s  state  is  somewhat  better  than  that  given  in  Table  1. 

21 
Goscinski   has  shown  that  Miller's  equation  can  be  derived  from  an 


h2 


expression  for  the  bracketing  function  ^   using  a  truncated  resolvent 
T  (p)  as  is  given  in  (l+.l)  and  that  the  points  £=£)t,  correspond  to 
its  solutions.   By  making  this  truncation,  however,  one  reduces  the 
contribution  from  the  continuum,  a  problem  that  is  eliminated  in  using 
the  Lbwdin  projection. 

A  study  was  made  to  determine  the  behavior  of  the  lower  bound 
curves  when  the  scaling  parameter  77  i-S  fixed  and  £»  is  permitted  to 
vary  over  a  wide  range  of  values.   For  both  He  and  Li  ,  g  was  taken  to 
be  4  and  NINT  to  be  hO.      In  Chapter  IV  we  observed  that  the  LBwdin  pro- 
jection greatly  changes  the  character  of  the  lower  bound  curves  from  that 
shown  in  Fig.  2  for  example.  We  found  that,  using  this  projection,  a 
minimum  appears  in  one  of  the  lower  bound  curves  at  each  of  the  points 
where  ^  equals  an  eigenvalue  of  ^   associated  with  an  eigenf unction 
in  the  reference  manifold.   This  behavior  is  shown  in  Figs.  3  a"^"^  ^« 
These  figures  demonstrate  the  effect  the  Lbwdin  projection  has  on  the 
eigenvalue  curves  of  £^  . 

It  is  interesting  to  examine  the  effect  that  changing  the  number 
of  internal  projection  functions  in  the  LBwdin  space  has  on  these  curves. 
Of  course,  using  fewer  functions  in  constructing  the  projection  t^ 
will  result  in  a  poorer  approximation  to  the  reaction  matrix  t^  . 
Keeping  ?7  =  1.2,  the  results  of  using  NINT  =  20  and  NINT  =  10  are 
shown  in  Figs.  5  and  6  respectively. 

We  observe  in  Fig.  5  that  the  lowest  bracketing  curve  Cii  <^oss 
not  recross  the  S  i  ~  £  line  after  the  minimum  at  £  =  E  ,  as  it  does 
in  Fig.  3  and  the  minimum  at   £)  =  E°  occurs  in  the  £)-,o  curve  and  not 
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in  the  ^  ^      curve  as  it  does  in  Fig.  3-   I^i  Fig.  o,  we  find  that  using 
only  10  functions  to  construct  the  LBwdin  projection  results  in  the 
minir.un  at  ^  =  "X  being  shifted  from  the   t  i -i  to  the  ^  ,_  curve. 
It  is  also  apparent  that  the   ^   ,  C  to  ^nd  C^ , ,  curves  vary  only 
slightly  from  the  unperturbed  eigenvalues  E^,  E  and  E,  respectively, 
over  the  region  under  consideration.   Only  the   C  -i -i  curve  is  appreciably 
affected  by  t^  . 

For  NINT  =  10,  we  find,  by  permitting  W  to  vary,  that  the  mini- 
mum for  2,  =  E   is  shifted  back  to  the  ^-,1  curve  for  Y)   =  1.35'   For 
^J  =   1.35  and  £  =    -2.06,  the  lower  bound  to  the  2  S  state  of  He  is 
-2.216  a.u.   These  results  are  shov.'n  in  Fig.  'J. 

We  conclude  this  study  of  the  lower  bound  curves  by  noting  that, 
if  the  LBwdin  projection  provides  an  adequate  approximation  to  the 
reaction  operator  t£  ,  the  i  th  lower  bound  curve  will  have  minimums  at 

each  point  ^=  E.,  where  E.  .  <E.  <  E.  and  will  cross  the  S  l    =  S 

J  1-1     J     1  1 

line  for  the  last  time  for  some  value  of  ^  not  an  eigenvalue  of  ^     . 
If  the  set  onto  which  the  projection  is  being  made  is  complete,  this 
latter  crossing  will  occur  at  the  eigenvalue  E.  of  CP  . 

This  behavior  is  observed  in  Fig.  3  where  the  &  •, -,  curve  has 
minima  for  ^  =  E^  and  ^=  E°  and  crosses  the  ^  ■!  =  ^  line  finally 
for  £  =  -2.1753.   From  (5.II)   £^3^  <  E°  <  E°  <  E  3^ ;  hence,  the  £   ^^ 
curve  should  have  minima  for  ^  =  E^,  E,  if  this  projection  adequately 
represents  the  reaction  operator.  We  see  in  Fig.  3  that  the  £  '   curve 
crosses  the  ^  •,  =  £  line  only  once,  at  the  minimum  £  l^  =   £  =  EO 


12   -    3' 
The  minimum  for  C  =  E,  appears  instead  in  the  P'   curve.   From  the 

^  l^       rf  f  ]_3 

results  in  Table  1,  we  see  that  C.  -i  q  P^^ovides  a  much  poorer  lower  bound 
to  the  3-"S  state  than  iS  1 -i  does  to  the  2"'S  state  of  He. 
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A  similar  analysis  can  be  made  of  the  Li  curves  shown  in  Fig.  k 
and  the  He  curves  shown  in  Figs.  5»  6>  ^^id  J, 

Although  the  choice  of  functions  of  the  type  (5 '9)  used  to  con- 
struct the  Lbwdin  projection  is  perhaps  not  the  best  possible  selection, 
tho  bracketing  curves  obtained  using  this  choice  illustrate  the  general 
behavior  we  may  expect  when  using  this  projection  and  the  material 
discussed  in  Chapter  IV. 


APPENDIX 

This  is  an  outline  of  the  procedures  used  in  the  calculations 
of  Chapter  V. 

The  calculation  of  the  LHwdin  projection  (U.22)  and  the  matrix 
^  ^  (i4-.13)  require  the  evaluation  of  the  quantities  ^t    j.,  i=l,2,  ..., 
I4O,  and  the  ^   matrix  whose  matrix  elements  are  given  by 

<Rj.iY-'/R^^  +  i  (€-e-)  <j,i9r>  <5>:/^;> . 

where  g  =   h. 

The  Hamiltonian  S^   and  the  inverse  of  the  perturbation  V  , 

'  o 

defined  by  (5.3)>  ™ay  be  expressed  in  terms  of  the  corrdinates  S,  t 
and  u: 

V'  =   2<.    »  (A.2) 
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From  (a. 3)  and  (5«9)  we  see  that  V   operating  on  the  functions 
j.  gives  a  linear  combination  of  terms  having  the  form 

where  n  =  -1,  0,  k  ^  0,  ^g  ^  0  and  m  >^  -1. 

From  this  we  can  see  that  the  matrix  elements  of  (A.l)  involving 
only  functions  of  the  LBwdin  space  are  simply  linear  combinations  of 
integrals  of  the  type 


Uu^)^  i'w'je"''stisjtrdujcs'-fr  t"j)t . 


(A.5) 


where  n  =  -1,  0,  1,  k  »  0,  ^  *  0  and  m  ^  0.  These  integrals  may  all  be 
evaluated  in  closed  form:  -^ 

^'  (A.7) 

_J I     ■ 

(^-IM^^  )-  -27r  "jTzp^sTT;^/  ^  Urr^f^^    Z^  I     v       (a.s) 
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The  integrals  involving  functions  of  the  reference  manifold  and 
the  LBwdin  space  ^(p.  [j,  /  can  be  expressed  as  a  linear  combination  of 
integrals  of  the  type 

"*  ^      '^  (A.9) 

where  i  ^  0,  k  ->  0  and  m  >  0.   The  solution  to  this  type  of  integral  is 
also  known  in  closed  form  and  is  given  by 

r~  D 


y^       ~     *"     Jle^   f.'     J  + 


^2-^.'  J  (A. 10) 
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The  number  of  integrals  of  this  type  and  of  the  type  given  by 
(a. 6),  (a. 7)  and  (A. 8)  was  quite  large,  so  a  program  was  written  for 
evaluating  them  on  the  IBM  7O9  computor. 

Once  the  integrals  were  solved,  the  matrix  elements  ^i^n  were 

calculated  and  the  ^  matrix  constructed.   The  inverse  of  this  matrix, 

-1  26 

/^  ,  was  obtained  using  the  method  of  successive  partitioning.    It 

was  then  a  rather  easy  task  to  determine  the  matrix  elements  of  ^-i  > 

which  are  given  by  the  formula 

Try  **t 

(A. 11) 

The  secular  equation  (4.30)  was  solved  using  the  subroutine 

JACFUL,  that  generates  the  eigenvalues  of  a  real  symmetric  matrix  by 

27 
the  Jacobi  method.  ' 
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